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Abstract 

The detailed mathematical model of heat and mass transfer of steel ingot of eurvi- 
linear continuous casting machine is proposed. The process of heat and mass transfer 
is described by nonlinear partial differential equations of parabolic type. Position 
of phase boundary is determined by Stefan conditions. The temperature of cooling 
water in mould channel is described by a special balance equation. Boundary con- 
ditions of secondary cooling zone include radiant and convective components of heat 
exchange and account for the complex mechanism of heat-conducting due to airmist 
cooling using compressed air and water. Convective heat-transfer coefficient of sec- 
ondary cooling zone is unknown and considered as distributed parameter. To solve 
this problem the algorithm of initial adjustment of parameter and the algorithm of 
operative adjustment are developed. 

1 Introduction 

Improved computing significantly increased role of mathematical modeling in research of 
thermo-physical processes. This, in turn, imposes stricter requirements towards accuracy 
and efficiency of mathematical models. 

It is well known that successful modeling mostly depends on the right choice of a 
model, which is directly affected by reliability of thermo-physical parameters used. Fre- 
quently, empirical data alone can not provide sufficient information about one-valuedness 
conditions. 

Therefore recently the big attention is given to the solution of inverse problems of heat 
conduction, in which it is necessary to define thermophysical properties of an object on 
available (frequently rather limited) information about temperature field. In particular 
thus it is possible to identify boundary conditions. There are difficulties in choice of some 
parameters of process for development of mathematical models of technological processes. 



While modeling process for specific industrial conditions it is necessary to determine 
some thermal or physical parameters each time, in particular convective heat-transfer 
coefficient (CHTC) on a surface of an ingot in the secondary cooling zone which depends 
on many factors. It is connected by that the convective heat transfer coefficient value is 
influenced with set of various factors. Besides, CHTC value can vary strongly enough in 
a time and on space coordinates. Thus, there is a problem of identification of the CHTC 
as distributed parameter. 

In the given work algorithms of initial adjustment of parameter when at the disposal 
of there is enough plenty of points in which the temperature on a surface of an ingot is 
measured, and operative adjustment when the temperature is measured only in one point 
on a surface are considered. 



2 Statement of problem 

The thermal field of the moving steel ingot and mold wall in the system of coordinates 
attached to motionless construction of CCM is considered [1]. In fig. 1 the diagram of 
CCM is introduced. 
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Figure 1: 



The heat conduction in the steel ingot in the mold area is described by nonstationary, 
nonlinear heat and mass transfer equation: 
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where v{t) - withdrawal rate, 21 - ingot thickness, Z - height of ingot in the mould, 
T(r, x, z) - metal temperature, c(T, x, z) - metal specific heat, p{T, x, z) - density, A(T, x, z) 
- thermal conduction, 6 — effective thickness of air gap between ingot and the mould wall, 
Xgz - thermal conduction coefficient of gap gas mixture, Tj^^^ - surface temperature of the 
ingot, T\^^i^^ - surface temperature of mold wall, (T„ - the resulted radiation coefficient. 

Conditions of equality of temperatures and Stefan conditions, and also boundary and 
initial conditions for the phase boundary are set: 
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where ^ - the phase boundary function of two variables x = ^{t, z), fi - crystallization 
latent heat. T/^,. crystallization temperature (average of the interval "liquidus - solidus"), 
n - normal to the boundary of phases. 
Heat equation for mould walls: 
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Boundary conditions for mould walls represent the character of heat exchange on each 
sight of wall: 
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where d - mold wall thickness, zq - mold wall altitude over meniscus level, ai - heat 
transfer coefficient from the mould wall to cooling water, T^ateriT, z) - cooling water 
temperature in the mold channel, q;2,3,4 ~ heat transfer coefficients from other mould wall 
to environment, Tos.2,3,4: - environment temperature, C„ - the resulted radiation coefficient. 

The following balance equation describes distribution of cooling water temperature in 
the mold channel: 



C - S ■ Vyjater """^^ = PjOll {Twater(,r, z) - T\^^^) - PsaE {Tyjaterir, z) - Te) , (6) 

where c - volume heat capacity of water, S - the cross-section area of the mold channel, 
Vwater - water velocity, Pj - perimeter of the interior mold wall, Pe - perimeter of the 
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external mold wall, - heat transfer coefficient from cooling water to the external mould 
wall, Te - external mould wall temperature. 

The cooling water temperature on the entry in the mould channel is known: 



and it's initial distribution in the mold channel: 
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The following equation describes heat and mass transfer on the curvilinear sections of 
CCM: 
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where 9m - angular velocity of ingot driving on the m-th curvilinear section. 
The conditions for unknown boundary on the curvilinear sections are 
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where ^i{ip) and 6(¥') ~ phase boundaries (interfaces). 

The boundary conditions of the secondary cooling zone include radiant and convective 
components of heat exchange and account for the complex mechanism of heat-conducting 
due to air-mist cooling using compressed air and water. The boundary conditions on the 
curvilinear sections are 
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where aj{G,n{T),(p), aE{Gm{T), ip) - convective heat transfer coefficients, Ci^,Ce^ - 
the resulted radiation coefficients, Ti^,Te^ — environment temperatures, Gm{T) — water 
discharge on the m-th section. 

The following equation describes the heat and mass transfer on rectilinear sections of 
CCM (analogously (1)): 
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When the liquid phase passes the straightening point on the rectilinear section of the 
secondary cooling zone, the conditions for the unknown phase boundary are set: 
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The boundary conditions for the rectilinear section: 
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We assume, that the thermal stream of the end of the rectilinear site is equal to zero: 
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The initial conditions for all temperature field (on the rectilinear and curvilinear sec- 
tions): 
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It is required to define the convective heat transfer coefficients a/(Grn(''"), <^), and 
«e(Gto,(t), (/:') using the available information about ingot temperature. 
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This is a boundary inverse problem and it is ill-posed in classical sense. Well-posedness 
in classical sense (or Hadamard wcll-poscdness) means performance of three conditions: an 
existence of a solution, its uniqueness and stability (input data continuous dependence). 
In our case the third condition is not satisfied. This is easily to verify using for the solution 
this problem the method of direct reversion [2] . Therefore other approaches are necessary 
to solve this problem. 



3 CHTC identification by least squares method 

Consider an ingot in first cooling section of secondary cooling zone. We have ingot surface 
temperature measurements in some points. So we have to solve the Dirichlet problem for 
interior heat exchange. The finite-difference method was used to approximate the solution 
of this problem. The convective heat-transfer coefficient (CHTC) has special distribution 
along the surface of the ingot. Parabolic function with a sufficient degree of accuracy 
approximates distribution of CHTC on the part of surface that is exposed to water-air 
spraying from one nozzle. This parabola has maximal value in the point that corresponds 
to nozzle coordinate. CHTC is considered as constant on the parts of the surface not 
subjected to the forced cooling (fig. 2). 




In one cooling section the same type spray nozzles are installed. They give an identical 
water-air spray. Hence the CHTC is the same parabola shifted along the abscissa axis (fig. 
2). 

All sites under spray nozzles can be reduced to the coordinate origin so that the peak 
of each parabola should be over the coordinate origin. Hence, it is necessary to define only 
two parameters - ap and Uc- So, a((/?) is given by 
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a{ip) = ac — ^(p'^ + a 



(18) 



Consider the parts of the section, on which a((p) = ac = const. Let K be the ensemble 
of points ipi, in which CHTC is equal to constant. Let B be the ensemble of other points. 
The finite-difference approximation of boundary condition (11) is 
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Ko ^''' ^^^ + ^^-'° = a,{Tj, - r,o) + Cj, {Tl - I^o), (19) 

where q - step of finite-difference grid by radius ri [3]. 

It follows that the discrepancy of heat flows on the boundary is: 
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Then we find a value ac, such that the sum of squares of discrepancies is minimum, 
i.e. the follow condition is satisfied 
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To the each point 99 j from we will put in conformity a point yi on the segment \—w^ w\ 
such that \yi\ is equal to the distance from the corresponding ipi to the coordinate of the 
nearest spray nozzle. Prom (18) and (19) we gain a discrepancy 
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On fig. 3 comparative results of calculations (1 - by the method of direct reversion, 
2 - by the least squares method) are presented. For steel grade st40, width of a slab 
is Im, 1 = 0,1m and v = l(m/minute). The decision obtained by the method of direct 
reversion is unstable and unsuitable for practical use. The second curve represents a spline 
approximation, which is gained as a result of the decision of a problem of identification by 
the least squares method. 
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Figure 3: 



Thus, we fined the spline approximation of the CHTC, which is distributed on the 
surface of the moving ingot. This approximation gives the minimum of mean-square 
deviation between measured surface temperature and calculated one according to the 
model as the result of solving of the direct problem. The CHTC for other sections of the 
secondary cooling zone is analogously defined. It should be noted that an advantage of 
the offered method is that the estimation error of the least squares method is negligibly 
small by relatively small number of abnormal measurements. It is very important in case 
of temperature measurement of a partially oxide scaled ingot surface. 
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4 Operative adjustment of convective heat transfer coeffi- 
cient (CHTC) 

CHTC obtained by initial adjustment varies under changes of various parameters of process 
(for example, ambient temperatures). Therefore, it is necessary to provide its operative 
adaptation during work CCM. The fine-tuning of parameters should be carried out in real 
time. But during usual work of CCM the information on a thermal condition of an ingot 
is limited to temperature indications in small number of points of the surface of an ingot. 
Such algorithms can be based on the stochastic approximation method [4]. 

The temperature on the ingot surface is measured in every equal small time intervals. 
Let us denote the measuring temperature data T* . The computer models the casting 
process using the presented mathematical model. The under model calculated temperature 
in the corresponding point we denote by Tj . It is necessary to correct the model parameters 
using information about deviations between measured and calculated temperature data to 
reduce these deviations to minimum. The difficulty of the decision of the given problem 
is that temperature measurements are deformed by a random telemetry error. 

Operative fine-tuning consists in refinement of the constant value ac, which defines 
the distribution of the convective heat transfer coefficient obtained by the solving of the 
problem of the initial adjustment of parameters. 

For using the algorithm of stochastic approximation it is necessary, that the random 
error of temperature indications would have the zero average and the finite variance. 

The algorithm of parameter adjustment is 

aj+i = aj - kj{T* - Tj), (20) 

where aj - j-th approximate value of ac, kj - special sequence of numbers, which satisfies 
to the following conditions: 

oo oo 

lim kj = 0, kj = oo, (kj)^ < oo. (21) 
For example the following elementary sequence satisfies to such conditions 

b + j 

where a, b € R, a > 0. Selecting numbers a and b, and also other sequences satisfying to 
the conditions (21), it is possible to change speed of convergence of algorithm. In [3], for 
example, it is recommended to keep kj as constant while the sign of discrepancy Tj — Tj 
not vary, and change then kj so that to satisfy to above mentioned restrictions. 

Truncation condition of the parameter fine tuning algorithm work is occurrence of m 
last received approximations an+i, an+2, , On+m in a vicinity of an serves: 

\an-an+i\<e, \/i = l,,m. 
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If the condition is executed, assume Oc is equal For check we use values CHTC 
which have been picked up experimentally at the decision of a direct problem of modeling 
of thermal field CCM [1]. 

5 Examples of realization of the stochastic approximation 
method 

Numerical modeling allows establishing the basic features of trajectories of parameter fine- 
tuning process. On fig. 4 trajectories of parameter fine-tuning, characterizing a deviation 
of the distributed parameter from true value, for the algorithm using sequence 

kj = j = 1, 2, 3, ... 

are presented at various values of factor a. When a < 1 very slow convergence is observed. 
In this case the time of parameter tuning is inadmissible big. 




t^emetry frame number 

Figure 4: 

If to choose a = 1 the value of the parameter is in enough small vicinity of true 
value approximately after 200th iteration. At a = 2 the trajectory of parameter fine- 
tuning reflects oscillations with damped amplitude and frequency and not later than for 
200 iterations the parameter is adjusted. At increase a > 2 the amplitude of oscillations 
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grows. In this case also oscillations with damped amplitude and frequency are observed, 
but for fine-tuning it is required considerably more iterations. 

From here we conclude, that for the chosen sequence the best values of the factor a is 
a number from interval 1 < a < 2. 

We investigate now influence of value b on speed of the algorithm's convergence. On 
fi-g. 5 trajectories of parameter fine-tuning are shown for various values b. Values b less 
than zero lead to that fine-tuning go in a "wrong" direction while the denominator is 
negative and at i = —b the denominator is equal to zero. Increase b leads to decrease of a 
velocity of convergence of algorithm. The same results have been obtained for sequences, 
which will be described further. Therefore further parameter b everywhere will be chosen 
equal to zero. 
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Figure 5: 

The following sequence also satisfies to conditions (21) 

'~ni ^^+i-\i + i, (i7-T,)(77+i-r,+i)<o • ^''^ 

Results of this algorithm work are presented on fig. 6. In this case factor a needs to 
be chosen within 1 < a < 3. Values out of this range give smaller speed of algorithm 
convergence. 

Consider another sequence, which also satisfies to conditions (21) 
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' - n, ' " 1 + 1, (77 - T,){T;^, - T,+i) < • ^''^ 

It has slower convergence than the previous two sequences. Results of calculations with 
use of this sequence are presented on fig. 7. Factor a can be chosen within 0.5 < a < 2. 
And, if 1.2 < a < 1.5, than obtained approximations differ from the true value no more 
than on 6 % after 20 iterations already. 




Figure 6: 

In the conclusion also it is necessary to notice, that the advantage of stochastic ap- 
proximation algorithm is its successful work for enough wide interval of initial values of 
the distributed parameter. 
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